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In 1929, Hermann Weyl derived m the massless solu¬ 
tions from the Dirac equation - the relativistic wave equa¬ 
tion for electrons. Neutrinos were thought, for decades, 
to be Weyl fermions until the discovery of the neutrino 
mass. Moreover, it has been suggested that low en¬ 
ergy excitations in condensed matterllHIl can be the 
solutions to the Weyl Hamiltonian. Recently, photons 
have also been proposed to emerge as Weyl particles in¬ 
side photonic crystals 10 • In all cases, two linear dis¬ 
persion bands in the three-dimensional (3D) momentum 
space intersect at a single degenerate point - the Weyl 
point. Remarkably, these Weyl points are monopoles of 
Berry flux with topological charges defined by the Chern 
numbers||2l|3|. These topological invariants enable mate¬ 
rials containing Weyl points to exhibit a wide variety of 
novel phenomena including surface Fermi arcsfl®, chi¬ 
ral anomalv ifITi . negative magnetoresistance[T2L nonlocal 
transportfO), quantum anomalous Hall effectfljl, uncon¬ 
ventional superconductivity! 151 and others llISlIlTl . Nev¬ 
ertheless, Weyl points are yet to be experimentally ob¬ 
served in nature. In this work, we report on precisely such 
an observation in an inversion-breaking 3D double-gyroid 
photonic crystal without breaking time-reversal symme¬ 
try. 

Weyl points are sources of quantized Berry flux of +2;7r in 
the momentum space. Their charges can be defined by the 
corresponding Chern numbers of +1, as shown in Fig. [T^. 
So, Weyl points robustly appear in pairs and can only be re¬ 
moved through pair annihilation. Since the Berry curvature 
is strictly zero under PT symmetry, — the product of par¬ 
ity (P, inversion) and time-reversal symmetry (T), isolated 
Weyl points only exist when at least one of F or T is broken. 
In Ref. (9), frequency-isolated Weyl points were predicted in 
FT-breaking DG photonic crystals. We chose to break P in¬ 
stead of r, in the experiment, to avoid using magnetic mate¬ 
rials and applying static magnetic fields. This also allows our 
approach to be directly extended to photonic crystals at opti¬ 
cal wavelengths. This F-breaking DG is shown in its body- 
centered-cubic (bcc) unit cell in Fig. [^. At the presence of 
r, there must exist even pairs of Weyl points. The two pairs 
of Weyl points illustrated in the Brillouin zone (BZ), in Fig. 
[T]:, are thus the minimum number of Weyl points possible. 
The bandstructure plotted in Fig. ^ shows two linear band- 
crossings along T - N and T - H. The other two Weyl points 
have identical dispersions due to T. 

We work at the microwave frequencies around lOGHz for 
the accessible fabrication of 3D photonic crystal. The current 
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FIG. 1. Weyl points found in the bandstructure of a P-breaking gy- 
roid photonic crystal, a) The Weyl points are monopoles of Berry 
flux in the momentum space, b) The bcc cell of the DG with 
a P-breaking defect on the red gyroid, where ai = (-1,1,1)|, 
a 2 = (1,-1,1)| and aa = (1,1,-1)|. The (101) surfaces are high¬ 
lighted in green, c) The BZ of the DG photonic crystal in b). Four 
Weyl points were illustrated on the green (101) plane along Y-H and 
T-N, where // = (0,1,0)^ and // = (-0.5,0,0.5)^. d) The photon 
dispersions were plotted along the A - F - //. The Weyl points are 
the linear band-touchings between the 4th-and 5th bands. 


additive processes like 3D printing can hardly fulfill the ma¬ 
terial requirement of low-loss dielectrics with high-dielectric 
constants. In order to fabricate the two inter-penetrating gy- 
roids with subtractive processes, we open up each gyroid net¬ 
work by layers along the [101] direction with equal thickness 
of . The P-breaking defects are introduced in each layer of 
red gyroid in Fig. [^. 

We approximate each gyroid network by three sets of hole¬ 
drilling, as illustrated in Fig. in a unit cell of the body- 
centered-cubic bcc lattice. Similar methods of drilling and an¬ 
gled etching has been used in fabrication of 3D photonic crys¬ 
tals at microwave CSl and near infrared wavelengths (13. The 
three cylindrical air holes, of the blue gyroid, along x, y and £ 
go through (0, |, 0)(2, (0,0, \)a and (|, 0, 0)a respectively. All 
air holes have a diameter of 0.54a, where a is the cubic lat¬ 
tice constant. Gyroids approximated by this drilling approach 
have almost identical bandstructures as those defined by the 
level-set iso-surfaces in Ref. 0. 

The second (red) gyroid is the inversion counterpart of the 









2 



FIG. 2. Fabrication of gyroids by drilling periodic holes along x, 
y, z directions, a) Illustration in a bcc unit cell that a single gyroid 
structure can be approximated by drilling air-holes, b) The double- 
gyroid structure can be made by stacking layers along the [101] di¬ 
rection. The red and blue gyroids, being inversion counterpart, inter¬ 
penetrates each other. We breaking parity (P, inversion) by shrink¬ 
ing the vertical connections to thin cylinders for the red gyroid. c) 
Shown on the right, a total of 18 layers were stacked. A zoom-in 
view from top is shown on the left with a ruler (in centimeters) in the 
background. 


blue gyroid. To break P, we shrink one arm of the red gyroid 
at ^)ri to be a cylinder oriented along [101]. Shown in 

Fig. [^, the defect cylinder has a diameter of O.la and height 
of 0.2a. We separate the red gyroid by cutting (101) planes at 
the centers of the defect cylinders. The blue gyroid, without 
introduced defects, is separated at the corresponding position 
of (-|, |, -^)a. Each layer is ^ thick, shown in Fig. m The 
separated layers of each gyroid are identical up to truncations 
of a 2 (the bcc lattice vector). One unit vertical period of DG 
consists of two layers from each gyroids, in which each layer 

is thick. 

V2 

The materials of choice are slabs of ceramic-filled plas¬ 
tics (C-STOCK AK from Cuming Microwave) of dielectric 
constant of 16 and loss tangent of 0.01. Each slab has a thick¬ 
ness of 9.5mm (= -^). Both width and length of the slabs 
are 304mm (~12 inches). The material hardness is adjusted 
between 80 to 90 on a Shore D gauge to be machined by car¬ 
bide tools without cracking. Computer-numerical-controlled 
machining was used to drill the slabs. Although all slabs of 
the same gyroid can be machined altogether, we processed no 
more than two layers at the same time to leave a solid frame 
after tilted drilling. Each layer experienced around 700 drills 
along +45 degrees away from its normal ([101]) and about 
40 drills in y on both top and bottom sides of the slab. A 
frame of ~20mm were left blank on the four sides of the slabs 
for handling. We subsequently cut two sides of the frames of 
each gyroid for assembling. Illustrated in Fig. [^, the lay¬ 
ers of the two gyroids were offset by half of the slab thick¬ 


ness (4.25mm). A picture of the final assembled sample is 
shown in Fig. [^. It consists of a total of 18 layers, with 9 
layers from each gyroid stacked in alternating order. 

We performed angled resolved transmission measurements 
on the photonic crystal sample to probe the dispersions of the 
bulk states. The schematic of the experimental setup is shown 
in Fig. [^. A pair of lens antennas were placed on the two sides 
of the sample. Transmission amplitudes (S21) were recorded 
by the network analyzer. The polarizations of both antennas 
was horizontal and their half-power beam widths were 9°. 
The collimated beam impinged on the sample (101) surface, 
in which the incident angle is varied by rotating the sample 
around the vertical axis. As illustrated in Fig. the com¬ 
ponent of the incident wavevectors [kosin(6)] parallel to the 
sample surface are conserved up to a reciprocal lattice vector 
of the sample lattice due to the discreet translational symme¬ 
try. All bulk states of the same wavevector projection (dashed 
line in Fig. [^) and frequency could be excited in the bulk and 
exist the sample in the same direction as the incident beam. 
The exiting signal was collected by the antenna at the oppo¬ 
site side. 

Illustrated in Fig. |^, the Weyl points in the system locate 
along G-H and G-N directions about half way between the 
BZ center (F) and the zone boundaries. To access them from 
free space, the incident angle 6 would have to be large than 60° 
and the effective cross-section of the sample would to small 
for enough signal to go through. To overcome this problem, 
we placed a pair of angled prisms (12.4°) in contact with the 
opposite surfaces of the sample as shown in Fig. [^. The 
prisms are made of the same material as the sample (e = 16). 

We mapped out all bulk states, projected along [101], by 
varying 6 and a shown in Fig. and b. Six transmission data 
of representative directions (a) are plotted in Fig. |^. The fig¬ 
ure insets are projected dispersions of the bulk bandstructures 
along the [101] direction. The transmission data stops on the 
right slanted boundary, which corresponds to the maximum 
rotation angle in 0. Close to this boundary, the transmission 
intensity is low due to the smaller effective cross-section of 
the samples at large angles. 

When a = 90°, the beam scan through the upper Weyl point 
along G-H represented by a magenta sphere. The data clearly 
shows a linear point touching at 11.3 GHz in frequency and ^ 
in wavevector. As a deviated from 90° to 105° and 120°, the 
point touching opens a gap as expected for a point-degeneracy. 
The other Weyl point, in cyan, on the right of the G - axis 
was studied by orient a to be 0°, 15° and 30°. Although the 
transmission intensity of the upper bulk bands is not promi¬ 
nent, the Weyl point dispersions can still be inferred from the 
curvatures of the lower bulk bands. In principle, all bulk bands 
of frequency and momentum matching the incident beam can 
be coupled and transmitted. However, the coupling and trans¬ 
mission efficiency depend on the details of the Bloch mode 
polarization, field distribution, group velocities, radiation life¬ 
time and so on. The low transmission amplitude in some parts 
of the data at high frequency region is due to the mismatch 
between bulk and the incident waves in free space. The re¬ 
maining two Weyl points, at the opposite k locations, relate 
to the two measured Weyl points by T. They have the same 
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FIG. 3. Transmission measurement results, a) Schematic of the microwave transmission setup, b) The bulk states in the sample are exited by 
the incident wavevector parallel to the sample surface varied by angle 6. c) Transmission data as the sample is rotated along the [101] axis by 
angle a. Insets are calculated bandstructures projected along [101]; they are scaled to the same range and ratio as the measured data. 


projected band dispersions and same transmission pattern as 
the data shown in Fig. [^. All the transmission data, in Fig. 
[^, compares very well with the theoretical bandstructures in 
the insets except a few regions with low transmission power. 

Our experimental demonstration of Weyl points answers 
the long standing search for a natural realization of the Weyl 
equations. These photonic Weyl points pave the way to topo¬ 
logical photonics 120^2^ in 3D, where 3D Dirac points 1241 
and various gapped topological phases 1^ can be discovered. 
Similar approaches can be readily adopted to observe Weyl 
points at optical frequencies using 3D nanofabrication ||26l 
[771 . We also anticipate that the Weyl points should be ob¬ 
servable in other systems [[28l, notably the topological semi¬ 
metals |[29l[30l in condensed matter. 
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